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ABSTRACT 

GRB 051 103 is considered to be a candidate soft gamma repeater (SGR) extragalactic giant 
magnetar flare by virtue of its proximity on the sky to M81/M82, as well as its time history, 
localization, and energy spectrum. We have derived a refined interplanetary network localiza- 
tion for this burst which reduces the size of the error box by over a factor of two. We examine 
its time history for evidence of a periodic component, which would be one signature of an 
SGR giant flare, and conclude that this component is neither detected nor detectable under 
reasonable assumptions. We analyze the time-resolved energy spectra of this event with im- 
proved time- and energy resolution, and conclude that although the spectrum is very hard, its 
temporal evolution at late times cannot be determined, which further complicates the giant 
flare association. We also present new optical observations reaching limiting magnitudes of 
R > 24.5, about 4 magnitudes deeper than previously reported. In tandem with serendipi- 
tous observations of M8 1 taken immediately before and one month after the burst, these place 
strong constraints on any rapidly variable sources in the region of the refined error ellipse 
proximate to M81. We do not find any convincing afterglow candidates from either back- 
ground galaxies or sources in M81, although within the refined error region we do locate two 
UV bright star forming regions which may host SGRs. A supernova remnant (SNR) within 
the error ellipse could provide further support for an SGR giant flare association, but we were 
unable to identify any SNR within the error ellipse. These data still do not allow strong con- 
straints on the nature of the GRB 051 103 progenitor, and suggest that candidate extragalactic 
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1 INTRODUCTION 

Giant flares are the most spectacular manifestations of soft gamma 
repeaters (SGRs). Their time histories are characterized by a very 
rapid (< 1 ms) rise to an intense peak lasting several hundred mil- 
liseconds, followed by a weaker, oscillatory phase which exhibits 
the period of the neutron star. The spectrum of the peak is very hard 
and extends to MeV energies. The most energetic giant flare to date 
is that of 2004 December 27 from SGR 1806-20, with an isotropic 
energy of well over 10 46 erg. (Hurley et al. 2005, Palmer et al. 2005, 
Mereghetti et al. 2005, Terasawa et al. 2005, Frederiks et al. 2007b). 
The observation of this event raised the question once more of the 
existence of extragalactic giant magnetar flares, and their relation 
to the short cosmic gamma-ray bursts (GRBs). Viewed from a large 
distance, only the initial peak of a giant flare would be detectable, 
and it would resemble a several hundred millisecond long, hard 
spectrum GRB. The energetics of giant flares make it a virtual cer- 
tainty that such events can be detected in distant galaxies, but rec- 
ognizing them and demonstrating their origin beyond a reasonable 
doubt remain difficult tasks. Predictions of the rates of extragalactic 
giant flares, and the percentage of short cosmic bursts which might 
actually be giant flares, vary widely (from a few to ~ 15% - Lazzati 
et al. 2005, Nakar et al. 2006, Popov and Stern 2006, Tanvir et al. 
2005), in part because of their unknown number-intensity relation 
(no SGR has yet been observed to emit more than one giant flare, 
and their distances are uncertain in most cases). However, they gen- 
erally agree that the rate is small enough that the majority of short 
bursts are indeed not due to SGR giant flares. On the other hand, 
the rate is not expected to be zero, so it is important to examine all 
possible cases exhaustively. In this paper, we revisit GRB 051103, 
a short burst whose location, time history, and energy spectrum are 
consistent with an origin as a giant flare in M8 1 . We make use of the 
full interplanetary network (IPN) data set to obtain a refined local- 
ization (an error ellipse). We utilize the capability of the RHESSI 
spacecraft to obtain time-resolved energy spectra with good energy 
resolution, at a time resolution which is limited only by counting 
statistics, and we take advantage of the fact that three instruments 
recorded energy spectra with good statistics to obtain joint fits. Our 
analysis follows that of Frederiks et al. (2007a), which was based 
on the Konus-Wind data. 

We also present new, much deeper optical data than previously 
reported for the section of the refined error ellipse closest to M81, 
taken 3 days after GRB 05110 3 (and approximately 16 hours af- 
ter the iGolenetskii et al.ll2005L GCN notice). We use these data to 
search for possible optical counterparts of this short burst (SGRB), 
and discuss the implications of its non-detection for its progenitor 
and putative association with M8 1 . Throughout this paper, we uti- 
lize the distance modulus of M81 determined by H ST observations 
of Cepheids, 27.8 (3.6 Mpc, iFreedman et"ail 19941) . 



2 IPN OBSERVATIONS AND LOCALIZATION 

GRB 051 103 occurred at 09:25:42 UT at Earth, and was observed 
by HETE-FREGATE (Atteia et al. 2003), RHESSI (The Ramaty 
High Energy Solar Spectroscopic Imager - Smith et al. 2002), 
Suzaku-WAM (Yamaoka et al. 2009), and Swift-BAT (Gehrels et 
al. 2004) in low-Earth orbit; the burst was outside the coded fields 
of view of Swift-BAT and the HETE-II imaging instruments, and 
was therefore not localized by them. It was also observed by INTE- 
GRAL SPI-ACS (Rau et al. 2005) at 0.5 light-seconds from Earth, 
Konus-Wind (Aptekar et al. 1995) at 4.5 light-seconds from Earth, 



and by Mars Odyssey (HEND and GRS experiments, Hurley et al. 
2006) at 232 light-seconds from Earth. A preliminary IPN error 
box was announced in a GCN Circular, where it was pointed out 
that this event had the largest peak flux ever observed by Konus- 
Wind for a short burst (Golenetskii et al. 2005). Optical follow-up 
observations were reported by Lipunov et al. (2005a,b), Klose et al. 
(2005), and Ofek et al. (2005, 2006), and a radio observation was 
reported by Cameron and Frail (2005). All yielded negative results. 
A MILAGRO GeV/TeV observation similarly produced only upper 
limits (Parkinson et al. 2005). 

The observations of each statistically independent spacecraft 
pair can be analyzed to produce an annulus of location, and the an- 
nuli can be combined to yield an error ellipse using the method de- 
scribed in Hurley et al. (2000). In this case, we have combined the 
Konus-HETE, Konus-RHESSI, Konus-INTEGRAL, Konus-Swift, 
and Konus-Odyssey annuli. Under these conditions (several rela- 
tively short baselines and one relatively long one), the error ellipse 
has a long major axis, corresponding to the annuli derived from the 
closer spacecraft pairs, and a short minor axis, corresponding to the 
annulus from the distant spacecraft pair. We obtain a 3 a error el- 
lipse centered at a(2000)=147.933° , <5(2000)=+68.589° , with ma- 
jor and minor axes 137'and 0.96'respectively, and area 104 square 
arcminutes. The chi-square for the error ellipse center is 0.9 for 3 
degrees of freedom (5 annuli minus 2 fitted coordinates). The area 
of the initial error box was 240 square arcminutefl The initial error 
box and the final error ellipse are shown in figure 1. 



3 TIME HISTORY 

The RHESSI time history of GRB 05 1 103 is shown in the top panel 
of figure 2. A distinctive signature of all three previously observed 
giant SGR flares within our Galaxy and the LMC to date is the 
periodic extended component following the initial short-duration 
peak. Among these three events, the periods of this extended tail 
have clustered around a narrow range of 5-8 seconds and also 
have a relatively narrow range of total isotropic energy releases of 
1 — 4 x 10 44 erg. This signal lasts for many minutes following 
the bursts but falls off rapidly after a few hundred seconds. While 
extended emission is frequently detected following cosmological 
short-hard bursts, such emission is not periodic. Therefore detec- 
tion of a periodic component of emission would be considered a 
strong confirmation of an SGR origin. 

None of the IPN light curves shows obvious evidence for ex- 
tended emission (pulsed or otherwise) following the burst. How- 
ever, it is conceivable that a marginally detected signal could be 
present within the noise. To search for such a component, we ac- 
quired Swift-BAT data for GRB 05 1 103 (binned at 64 ms) and used 
the Lomb (1976) periodogram to calculate the relative power in 
the signal following the burst at periods up to about 20 seconds. 
We created periodograms for all of the four BAT energy channels, 
which cover the energy range 15 to 350 keV (and for combinations 
of channel sums) and for various time ranges following the emis- 
sion (ranging from the first 60 s to the first 300 s.) To assess the 
significance of any peaks in the power spectrum, we performed a 
Monte Carlo analysis by repeatedly randomizing the order of the 
64 ms time bins for each data set over the range of interest and 



1 A typographical error in GCN 4197 incorrectly gave the area as 120 
square arcminutes 
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Figure 1. This shows in red the original error trapezium provided by the 
IPN, in green the 3<r refined error ellipse for the position of GRB 051103 
and in blue the fields of the region studied using KPNO. The asterisk in- 
dicates the center of the ellipse and the most likely arrival direction of 
the burst. Approximately 65 square arcminutes of the ellipse are contained 
within the old error box. These are imposed upon an image of the area sur- 
rounding M8I from the Digital Sky Survey. The possibility that this burst 
came from the inner disk of M8 1 is firmly ruled out. However, the bright- 
est GALEX UV knots noted by Ofek et al. (2006) are within the ellipse. 
Lipunov (2005b) noted the presence of two galaxies within the initial error 
box, PGC 2719634 and PGC 028505. The former galaxy lies at the 18% 
confidence contour of the ellipse, and remains a plausible host candidate, 
while the latter lies at the 0.03% contour, and is unlikely to be the host. 



measuring the rate of occurrence of independent peaks above var- 
ious power levels. We identified no peaks with greater than 98% 
significance in any channel or time range. 

This non-detection is expected. To assess the general de- 
tectability of periodic post-flare emission from extragalactic giant 
magnetar flares, we also acquired the Swift-BAT light curve of the 
the 2004 December 27th flare from SGR 1806-20. We then scaled 
the signal down by a factor of (D/Dsgr) 2 and added it to the light 
curve of GRB 051103 (both scrambled and unscrambled). No sig- 
nal is detected in the periodogram at the known periodicity of 7.56 
seconds at the distance of M8 1/82 (D = 3600 kpc). The maximum 
distance for detecting periodicity with our analysis greater than 3 
sigma is only about D = 220 kpc (if a distance of TJsgr = 14.5 
kpc to SGR 1806-20 is assumed, or to 130 kpc if the 8.7 kpc dis- 
tance of Bibby et al. 2008 is assumed), less than the distance even 
to M31. This limit may not be exact: both SGR 1806-20 and GRB 
051103 were detected off-axis by BAT and the comparative satel- 
lite sensitivity will depend on the specifics of the off-axis angle. 
For the giant flare from 1806-20, BAT was pointing 105°away, and 
slewed to 61°away starting around 38 s after the peak. For 051 103, 
BAT was pointing 122°from the source. However, as the expected 
signal from a December 27-like event at the distance of M81/M82 
would be only 0.01 sigma assuming similar sensitivities for the two 
events, we consider it extremely unlikely that any possible angle 
outside the BAT FOV would lead to a detection unless the periodic 
component were several orders of magnitude stronger than that ob- 
served in the three Galactic/LMC events to date. 



4 ENERGY SPECTRUM 

A key signature of the spectra of the three SGR giant flares ob- 
served to date is a very hard energy spectrum for the initial, sev- 
eral hundred millisecond long burst, and a dramatic spectral evolu- 
tion to a soft spectrum for the subsequent pulsating component. As 
these bursts were observed by various instruments, with different 
temporal resolutions, spectral resolutions, and energy ranges, and 
all of them were in some degree of saturation at the peak, a precise 
description of the spectra is impossible. Nevertheless, all of them 
can be characterized as very hard spectra at the peak, sometimes 
consistent with a very high temperature blackbody (e.g. Mazets et 
al. 1979, Fenimore et al. 1981, Hurley et al. 1999, Mazets et al. 
1999, Hurley et al. 2005, Frederiks et al. 2007b). Accordingly, we 
have analyzed the time-resolved energy spectra of GRB 051103. 
RHESSI, Konus, and Suzaku obtained energy spectra for GRB 
051103 over a wide energy range, with good statistics, although 
with different time resolutions. (Due to the off-axis arrival angles 
at Swift and HETE-II, the detector response matrices are not well 
known, and we have not used these data.) Because the finest time 
resolution can be obtained from the RHESSI data, we have ana- 
lyzed the RHESSI spectra both separately, to obtain the best time 
resolution, limited only by counting statistics, and combined with 
the Konus and Suzaku data, to obtain the best statistics, albeit at the 
cost of temporal resolution. 

RHESSI uses nine unshielded coaxial germanium detectors 
to observe a broad energy band (30 keV-17 MeV) with excellent 
energy (1-5 keV) and time resolution (1 binary /is) and moderate 
effective area (~150 cm 2 ). The data are recorded event-by-event, 
which provides great flexibility in choosing analysis intervals. 

To determine RHESSI's spectral response to GRB 05 1 103, we 
used the Monte Carlo package MGEANT (Sturner et al. 2000). We 
simulated monoenergetic photons in 192 logarithmic energy bins 
ranging from 30 keV-30 MeV generated along a 60°azimuthal arc 
at the 97°off-axis angle of GRB 05 1 103. We fit a polynomial back- 
ground and extracted the burst data in SolarSoftWare-IDL routing 
Because of radiation damage to some of the detectors, we used only 
data from rear segments 1, 4, 6, 7, and 8. Spectral fitting was con- 
ducted with ISIS v 1.4.9 (Houck 2000). In general, the full 30 keV- 
17 MeV energy band was employed, except when sufficient counts 
could not be accumulated at high energies. 

We fit the data with a Band (Band et al. 1993) function: 

N = J A{E/E ptv ) a exp(-E/E ) E < E break 

]^ B(E I Epiy) 13 E > Ebreak 

wilhE break = E (a-P) andB = A{ ^=^ ) a - fi exp(p-a). 
For f3 < —2 and a > —2, E pea k = Eo(2 + a) corresponds 
to the peak of the uF u spectrum. The normalization A has units 
photons/(cm 2 s keV), and E P i V is here taken to be 100 keV. For 
joint fits, the Band function parameters a, /3, and E pea k were tied 
for both instruments, but the normalizations were allowed to vary 
independently. 

For the RHESSI-only time-resolved fits, we identified time in- 
tervals with background-subtracted S/N of 20 in the 60 keV-3 MeV 
band. This yielded three intervals, to which a fourth tail interval 
of S/N =12 was added. For most intervals, the cutoff power law 
model (CPL), equivalent to the Band function below Ebreak, pro- 
vided the best fit. The time evolution of the parameters of the best- 
fit spectral model (a cutoff power-law) is presented in the lower 

2 http://www.lmsal.com/solarsoft/ 
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Figure 2. Time history of GRB 051103, and evolution of the spectrum. 
The top plot shows the dead-time corrected RHESSI light curve (60 keV- 
3 MeV) with 1 millisecond time resolution, starting at 09:25:42.184 UT. 
The background count rate is 0.55 counts/ms and has not been subtracted. 
The time history has an e-folding rise time of 1 .2 ± 0.04 ms, an e-folding 
decay time of 28.6 ± 0.6 ms, and a Tgo of 100 ± 4 ms. The middle and bot- 
tom plots show the evolution of the best-fit peak spectral energy and power- 
law index for the cutoff power-law model. The black points are RHESSI- 
only, while the gray points are joint fits between RHESSI and Konus-Wind. 
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Figure 3. Joint spectral fit of RHESSI (black) and Suzaku-WAM (gray) data 
for a one- second interval containing the burst. 



ence; this deficit appears to be a result of data lost due to deadtime 
during the intense peak of emission. 

Konus-Wind triggered on GRB 051103 and recorded 64 ms 
spectra in the 20 keV - 10 MeV range; we conducted joint fits 
between RHESSI and Konus-Wind for the 64 ms and 128 ms inter- 
vals analyzed in Frederiks et al. 2007a. These fits are presented in 
Figures 3 and 4, and the details are reported in table 1 . Good corre- 
spondence was obtained in the best-fit parameters between the two 
instruments, although a normalization offset was necessary. 

The spectrum of the 2004 December 27 giant flare from 
SGR1806-20 was measured by many different instruments, using 
many different methods (Hurley et al. 2005, Boggs et al. 2007, 
Palmer et al. 2005, Frederiks et al. 2007b). While they do not agree 
on the exact shape of the spectrum, none found evidence for the ex- 
istence of a high energy power law component in the Band model. 
Our RHESSI-only spectral fits of GRB 051103 are consistent with 
this, but in two of the joint fits this component is preferred (table 
1). A blackbody can be fit to the RHESSI data, but only over the 
800 keV - 5 MeV range; kT is approximately 550 keV for this fit, 
and the chi-square is about 1.5 per degree of freedom. 

Table 2 contains a comparison of the energetics of the giant 
flares observed to date. Because of the uncertainties in the SGR 
distances, as well as the different energy ranges, time resolutions, 
and characteristics of the various instruments which observed these 
events, comparisons between the values given are probably uncer- 
tain by at least a factor of 3. The total energy of GRB 051103 at 
the distance of M81, 7.5 x 10 46 erg, does not strain the possibilities 
of the magnetar model. However, PGC 2719634, whose distance is 
unknown, also remains a possible host (Lipunov et al. 2005b). 



panels of Figure 2. The initial spike of emission has a significantly 
higher peak energy than the decaying tail; however, the spectral in- 
dex of the power law appears to harden throughout the burst. The 
results are reported in table 1 . 

Suzaku-WAM did not trigger on GRB 051103, so the only 
data available are for a 1 -second spectrum containing the entire 
burst, in the 50 keV - 5 MeV energy range. The RHESSI- WAM 
joint fit is shown in Figure 3 and the fit results are reported in table 
1. The WAM fluence is a factor of ~5 lower than the RHESSI flu- 



5 OPTICAL OBSERVATIONS AND ANALYSIS 

Observations were obtained on 6 th November 2005 using the Mo- 
saic wide field optical imaging camera at the KPNO 4m tele- 
scope. These data reach a limiting magnitude of ~ 24.5 in the R 
band, which is co nsiderably deeper than the study completed by 
lOfek et all d2006h - The observations covered the majority of the 
original error region, and in particular that part nearest to the galaxy 
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Table 1. Best fit parameters for the joint fits of the RHESSI/Suzaku-WAM, RHESSI/Konus-Wind, and RHESSI-only data. 
Times are relative to TO = 09:25:42.190 UT in the RHESSI frame. Joint fit fluences are in the 20 keV-10 MeV band, while 
the RHESSI-only fluences are 30 keV-10 MeV. The instrument normalizations were free to float in the fit; the normalization 
of the second instrument relative to RHESSI is given. Errors are quoted at the 90% confidence level. 
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93.0/65 = 1.43 
30.9/32 = 0.96 
38.5/11 = 3.50 
7.5/11 =0.68 
3.7/4 = 0.93 
17.8/7 = 2.54 



Table 2. Approximate energies and peak luminosities of the SGR giant flares, and of GRB 051103. 



SGR Energy, erg Peak Luminosity, erg cm 2 s 1 Assumed distance, kpc 



0525-66^ 
1900+14^ 


1.2 x 10 44 


5 x 10 44 
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4.3 x 10 44 


2 x 10 46 


15 
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15 
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1.5 x 10 45 


1.2 x 10 47 


780 


GRB 051103 - 


7.5 x 10 46 


4.7 x 10 48 


3600 



a Mazets et al. 1979 

b Hurley et al. 1999; Tanaka et al. 2007 

c Hurley et al. 2005; Terasawa et al. 2005; Frederiks et al. 2007b 

d Mazets et al. 2008 




Figure 4. Joint spectral fit of RHESSI (black) and Konus-Wind (gray) data 
for interval 1. The overplot model is the best-fit Band function; the normal- 
ization between the datasets was allowed to vary in the fit. 



Figure 5. Joint spectral fit of RHESSI and Konus-Wind data for interval 2- 
3; symbols as in Figure |4] Two high -energy RHESSI points were negative 
and are omitted from the logarithmic counts plot. 



M81. The images were flat-fielded and sky subtracted using stan- p or comparison, pipeline-reduced images of the region from 

dard tasks within IrafQ. 

3 IRAF is distributed by the National Optical Astronomy Observatories, Astronomy, Inc., under cooperative agreement with the National Science 

which are operated by the Association of Universities for Research in Foundation. 



6 Hurley et al. 



CFHT LS were obtained via the Virtual Observatory ( Walt on et al.l 
2006). These formed part of the wide synoptic survey i n the R band, 
also to a limiting magnitude of ~25 dllbert et al.l2006h . Coinciden- 
tally, the region was imaged on I s * November 2005, 2 days prior to 
the burst, and re-imaged within 1 month after the burst. This pro- 
vided an ideal data set for comparison to the KPNO images as the 
timescale between the first images by CFHT and our images from 
KPNO is only 6 days, minimising any modulation in long-period 
variable stars in the disk/halo of M81. 

Figure 1 shows the previous error quadrilateral, the refined 
3a error ellipse and the fields covered by our KPNO observations, 
in relation to M81. Our observations were positioned to cover the 
original error quadrilateral but still cover 62% of the refined 3cr el- 
lipse and contain 76% of the total likelihood. It is important to note 
that our observations cover the region closest to M81, and therefore 
our search addresses the possible association of GRB 051 103 with 
M81. 

Initially, we searched the images for variability of afterglow 
counterparts, either at the distance of M81 or in the background, by 
visual inspection and no obvious afterglow candidate was found. 
The magnitudes of sources within these images were then studied 
using S EXTRACTOR within GAIA. They were all calibrated to the 
r band magnitudes of stars in the surrou nding region as published 
in the Sloan Digital Sky Survey (SDSS) dAdelman-McCarthv etal] 
2008). The r band filter used by CFHT matched the filter used in 
SDSS and the filter used by KPNO was a Cousins R band filter. Al- 
though this is partially taken into account in the calibration to r band 
magnitudes, there are some sources which have large colour differ- 
ences, for example very red sources. If a source appeared to differ 
in magnitude between the CFHT and KPNO images, th e colour 
correc tion was calculated using a formula developed by iLuptonl 
( 2005) and it was then determined if the magnitude difference was 
due to colour effects. If it was not due to colour effects, the source 
was investigated further. It is important to note that there may be 
a source within the field which was varying but has not been iden- 
tified due to this colour correction method. However, this method 
would only miss objects with a variability of ^0.3 magnitude (the 
average colour correction factor used). 

Although the magnitude of some stars differed between the 
images, these were found to be caused by other factors, for exam- 
ple, being near chip edges or large diffuse galaxies unidentified by 
SEXTRACTOR. One of the stars in the region studied has a vary- 
ing magnitude on the images studied and further investigation con- 
firmed it is likely a variable star. 

We checked extended sources to look for a conventional 
SGRB afterglow within a moderately distant host galaxy, with a 
limiting magnitude of ~ 23.3. If an extended source appeared to 
be varying due to a possible point source being superimposed on 
it, the colour correction was calculated and the object was studied 
in more depth by eye. This involved using the software to match 
seeing conditions and measure the size of the object, and then to 
check if there was an indication of a change in shape which might 
indicate a superimposed afterglow component. 

In addition to the photometry described above we also 
searched for afterglow candidates with PSF-mat ched image sub- 
traction, using a m odified version of the ISIS code l lAlard & Luptonl 
This method gives us a better chance of finding 
sources that are blended with other, brighter objects (i.e. bright host 
galaxies). After cosmic ray cleaning and resampling on a common 
pixel grid, we subtract the KPNO data with the CFHT data taken 
before and with the data taken after the burst. We found no credible 
afterglow candidates. 



Table 3. The observed fluence, in the energy band 15-150keV, of SGRBs 
with observed R band magnitudes at approximately 3 days. 



SGRB 


Fluence 
10 — 7 erg cm" 2 


R band Magnitude at 3 days 


051221A 


11.6±0.4 CO 


24.12±0.28 ( 2 ) 


051227 


2.3±0.3 < 3 ) 


25.49±0.09 C 4 ) 


060121 


2 6 7+5.3 (5) 


25±0.25 < 6 > 


060614 


217±4 ( 7 ) 


22.74±0.31 < 8 ) 


061006 


14.3±1.4 < 9 ' 


>23.96±0.12 ( 10 ) 


070707 


334+0.753 (11) 
u.jjt_ 316 


26.62±0.18 < 12 ) 


070714B 


7.2±0.9 < 13 ) 


<25.5 < 14 ) 


071227 


2.2±0.3 ( 15 ) 


>24.9 ( 16 ) 


080503 


20.0±1 (16 ) 


25.90±0.23 < 17 ) 



(0 ICummings et alj J2005I) < 2 > ISoderberg et alj fepod) < 3 ) iHullinger et al J 
12005J) < 4 ) iD'Avanzo et all d2009l) t s > an approximate fluence calculated 
using spectral para meters published bv[Golenetski i et alj d 2006h ^ b ased 
on observation s by Eevan et alj l2006ah ( 7 ) iBarthelmv et alj <2006l) ( 8 > 
iMangano et at] 120071) < 9 > iKrimm et alj 120061) CtO) an upper limit based 
on observations 2 days after the burst completed bvlD'Avanzo et al] J2009h 
C 11 ) an a pproximate fluenc e calcu lated using spectral parameters pub- 
lished by iGolenetskii et alj l2007t) t 12 ) Piranomont e et alj <2008t) ( 13 ) 
iBarbier et alj J2007h ^ 14 ^ a lower limit based on observations 4 day s af- 
ter the burst completed bv lPerlev et"ai] j2009h < 15 > |s a to et all j2007n (16) 
a 3(7 u pper limit published bylD'Avanzo et al] d2009l) ( 17 ' lUkwatta et alj 
120081) ( 18 > |PerlevetaTll2009l) 



The analysis of the images found no optical afterglow candi- 
date in the region studied 3 days after GRB 051 103. This can place 
constraints on the progenitor of GRB 051103 by considering the 
expected results for the potential progenitors. 



5.1 Progenitor option 1: a Short GRB (SGRB) 

The optical afterglows of various SGRBs have been studied and 
these data can be used to predict the range of afterglow properties 
of an SGRB of a particular gamma-ray fluence. There is evidence 
for a reasonable correlatio n, to first order, between gamma-ray 
fluence and afterg low flux dNvsewander. Fruchter. & Pe'erl 120081 ; 
iGehrels et al]|200a) . Using XSPEC, we created a model spectrum 
of GRB 051 103, using the RHESSI + KW joint fits in Table 1, and 
estimated the fluence of GRB 051103 in the energy band 15-150 
keV to be approximately 9.6ig 4 7 5 x 10 _7 erg cm -2 . 

It is possib le to compare GRB 05 1 103 to other SGRBs in the 
BAT catalogue dSakamoto etai]|2008l) using the approximate flu- 
ence, calculated for the energy band 15-150 keV, and the photon 
indices given in Table 1. GRB 051103 is isolated at the extreme 
bright, hard end of the SGRBs in the Swift distribution (c.f. Fig- 
ure 14 from Sak amoto et al. 2008). Similarly, in the study of short 
bursts by Mazets et al. (2004) over the much wider Konus energy 
range, of the 109 spectra which could be characterized by an E pea k, 
none exceeded 2.53 MeV. The peak energy of GRB 051103 is ap- 
proximately 3 MeV (table 1). Thus if GRB 051103 is an SGRB 
rather than an SGR giant flare, it is a fairly extreme case. 

We compared the fluence of this burst to other SGRBs ob- 
served by the Swift Satellite. Table 3 provides the data of SGRB 
with fluences in the band 15-150keV and late optical observations, 
obtained from the relevant GCNs, and measured optical afterglows 
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in the R band, approximately 3 days after each burst 0. For two of 
the bursts it was necessary to estimate the fiuence in the correct 
energy band using the same method as with GRB 051103. This is 
not a complete sample of SGRBs, as there are a number with a 
relatively low gamma-ray fiuence that were either not observed op- 
tically, were not observed for longer than a few hours, or did not 
have a detected optical afterglow. We chose this sample so we did 
not have to rely on the assumption that we can extrapolate the light 
curve to later epochs and because they are of a similar gamma-ray 
fiuence to GRB 051103. We compared the SGRBs in Table 3 to 
GRB 051103 and predict the optical afterglow would have an R 
band magnitude of ~ 24 as it is at the higher end of the fiuence 
distribution. This is within the limiting magnitude of the KPNO 
and CFHTLS images us ed, but would hav e been unobservable in 
the images obtained by lOfek et al.l J2006h . As no afterglow was 
observed, this rules out most typical SGRBs in the region of the 
error ellipse covered by our imaging. However, there are cases of 
SGRBs with extremely faint optical afterglows, for example GRB 
080503, which had a similar fiuence to GRB 051 103 an d an r band 
magnitude of 25.90±0.23 at 3 days jPerlev et alj|2009l) . So the ob- 
servations cannot rule out an unusually faint SGRB in this region 
similar to GRB 080503. Additionally, GRB 05 1 103 could be a clas- 
sical SGRB in the part of the error ellipse not studied in this paper. 

5.2 Progenitor option 2: an SGR giant flare in M81 

Conversely, GRB 051103 could be an SGR giant flare in 
M81 with similar ener gy to the giant flare from SGR 1 806-20 
dGolen etskii et al. 2005]) and a very faint optical afterglow (Eichler 
l2002t l iLevan et alj2008h . Using observations of the giant flare from 
SGR 1806-20, we can predict the apparent optical magnitude of an 
SGR in M81. The distance to SGR 1806-20 has proven difficult to 
determine; the dis tance mod ulus adopted by ma ny authors is 15.8 
dCorbeletal.1 19971) . although Bi bbv et al. (2008) recently obtained 
a revised distance modulus estimate of 14.7±0.35mag. Here we 
continue to use the larger distance modulus as this will provide an 
approximate upper limit on the absolute magnitude. The giant flare 
from SGR 1 806-20 had an obse rved radio afterglow and this has 
been used bv lWang et alj l l2005i) to make predictions of the appar- 
ent R band magnitude of the afterglow. Their analysis suggests that 
the giant flare would have had an apparent magnitude of ~ 22 at 
3 days, and hence an absolute magnitude of M w 6. Taking this 
as the absolute magnitude of any afterglow of GRB 051103 if it 
is an SGR giant flare, and using the distance modulus to M8 1 of 
27.8 jFreedman etaill 19941) . we conclude the afterglow would be 
expected to have an apparent magnitude of > 34. Despite the many 
uncertainties involved in this calculation, we can have some confi- 
dence that such an afterglow would not be detectable with the data 
available. For future reference, it is important to note that with more 
accurate positions and rapid follow up observations it may be pos- 
sible to observe the optical afterglows of extragalactic giant flares. 
For example, if there were a second potential giant flare in M81 we 
predict the optical afterglow would have a peak apparent K band 
magnitude of ~20 at 86s after the giant flare and would fall to ~26 
at 1 hour. This is observable with current and upcoming facilities, 
for example the European Extremely Large Telescope (E-ELT). 
However, these predictions are based upon the theoretical mod- 
els of SGR giant flares being similar to the blast wave model used 

It is important to n ote the classificat ion of some of these SGRBs are cur- 
rently being debated fchang et alj2009h 



to describe classical GRBs. IWang et alj J2005h use the blast wave 
model and radio observations of the giant flare from SGR 1806- 
20 to extrapolate the optical afterglow. S GRs have been obse rved 
during periods of act ivity using ROTSE-I {Akerlof et al . 2000) and 
Swift (for example ICummings et al.l l2009). and IR observations 
have been obtained fo r SGR 1900+14 4.1 days after outburst de- 
tecting no variability (Opp enheimer et aT] [l~998). These have pro- 
vided upper limits on the optical afterglows from the softer spec- 
trum, shorter, and weaker bursts seen during active phases of SGRs 
but it is important to note that there have been no reported rapid 
optical follow up observations of galactic SGR giant flares, which 
have a significantly higher fiuence and are spectrally harder than 
these bursts. Therefore, we are completely reliant on theoretical 
predictions and future observations may show discrepancies with 
these predictions. Indeed, our observations with a limiting magni- 
tude of 24.5, giving an absolute magnitude -3.3 assuming it is at 
a distance of 3.6 Mpc, constitute one of the deepest absolute mag- 
nitude searches for an afterglow from a possible SGR giant flare. 
This absolute magnitude is only exceeded by the search for an af- 
terglow from GRB 070201, which is a candidate SGR giant flare 
in M31, corresponding to an absolute magn itude of -7.4 obtained 
10.6 hours after the burst dOfek et alJ l2008). However, as we dis- 
cuss l ater, it is unlikely that both of thes e events were SGR giant 
flares (Chapman . Priddev, & Tanvij2009t) . 

From the GALEX UV imaging (Mar tin et al]l2005l) . there is 
evidence that the error ellipse does contain star forming regions in 
the outer di sk of M81. The t wo brightest UV sources are marked 
on Figure 6 dOfek et alj2006h . These young stellar regions in M8 1 
could host an SGR which could emit a giant flare. Similarly, these 
UV regions could be the locations of massive star clusters, and 
SGRs 1900 +14 and 1806-20 have been associated with massive 
star clusters ( Mirab el & Fuchsl 19991 ; Pyuria et alj200d) . However, if 
GRB 05 1 103 is an SGR giant flare in M81, we might also expect t o 
find a young (up to ~ 10 4 years old. iDuncan & Thompsonl dl992l) ) 
SNR in the nearby region, although this association is still being 
debated dGaensler et al j200ll . 120051) . When an SGR is formed, it is 
theore tically possible that it is give n a kick of up to 1000 km s _1 or 
more (Duncan &Thompsonlll992t) and therefore could have trav- 
eled a distance of >10pc from the SNR. However, this is only 
equivalent to an angular separati on of ~ 0.6arcsec at a distance 
of 3.6Mpc dFreedmanetaljl 19941) . Hence, an accompanying SNR 
would still be expected to fall within the error ellipse. Of the known 
SNR in M81 dMatonick & Fesenll 19971) . there are none within the 
error ellipse. 

M81 has bee n studied by the Chandra X-Ray Observatory 
dSwartz et alj2003l) and three X-ray sources are within the error el- 
lipse. However, they have not been identified in visible or radio ob- 
servations. Additionally, they have not been identified with known 
SNR, nearby stars, are not co-incide nt with HII starforming region s 
(the expected location of SGRs - IDuncan & Thompsonl dl9920 ) 
and ar e more likely to b e X-ray binary systems than unidentified 
SNR dSwartz et alj2003h . This survey had a limiting luminosity of 
3 x 10 36 erg s _1 , which means it would detect very young su- 
pernovae, as the oldest supernovae with detected X-ray afterglows 
had a luminosity of ~ 10 37 erg s _1 and an age of ~ 60 years 
dSoria & Pern3 [2008). Additionally, this survey would detect the 
X-ray luminous SNR as these have a luminosity of up to ~ 10 37 
erg s~ but would not detec t the X-ray faint SNRs which have a lu- 
minosity of ~ 10 34 erg s" 1 dlmmler & Kunt J2005I). SGRs are well 
known to be quiescent soft X-ray emitters and iMereghetti et al.l 
(2000) have measured the soft X-r ay flux of SGR 1806-20 to be 
approximately 10~ erg cm" " 2 s- MFrederiks et alj d2007ah deter- 
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mined that the Chandra Observatory would be unable to observe 
directly the persistent X-ray flux from an SGR like SGR 1806-20 
in M81. 

An alternative method of searching for SNR is to use Ha 
and OIII narrow band observations. The Isaac Newton Telescope 
(INT) has been use d to search for planetary nebulae in M81 by 
Magri nietaT] ( l200ll ) and they have found 171 potential candidates, 
some of which are in the nearby region of the refined error ellipse. 
Their criteria for differentiating between an SNR and a planetary 
nebula is that planetary nebulae cannot be spatially resolved and 
SNR are. A young SNR, as required for an SGR, could be misiden- 
tified as a planetary nebula by this criterion, since a one arcsecond 
region corresponds to a physical size of ~ 20pc. Young SNRs may 
well be significantly smaller than this, since an expansion veloc- 
ity of a few thousand km s" 1 over a magnetar lifetime of ~ 10 4 
years leads to sizes of 10 — 50 pc. Indeed, many S NRs in M82 
appea r (based on radio maps) to be fairly compact dFenech et al.l 
2008). However, the nearest is still ~23 arcsec from the error el- 
lipse, and we showed earlier that an SGR in M81 would only be 
able to travel ~0.6 arcsec from its birthplace. The Ha luminosity 
of SNR s in nearby disk gal axies tends to be greater than ~ 10 36 
erg s _1 dde Griis et al]2000|) and the work bv lMagrini et alj J200lh 
had a limiting Ha flux of less than 6 x 10" 17 erg cm" 2 s" 1 which 
corresponds to a limiting luminosity of ~ 10 35 erg s _1 . Therefore, 
we would expect their survey to find a candidate SNR. We used the 
recently published Ha an d OIII images from the INT Wide Field 
Camera Imaging Survey dMcMahon et aljfeooih . with a limiting 
luminosity of ~ 10 3 ° erg s - 1 at the distan c e of M 81 as these are 
the same images as used by iMagrini et al.l d200ll). and compared 
them with 21cm radio images from THINGS dWalter et alj|20o3) 
and Chandra X-ray source positions dSwartz et alj|2003h to search 
for previously unidentified SNRs within the error ellipse. There is 
a possible circular 21cm radio source coincident with a Chandra 
X-ray source of approximately the right flux for an SNR in M81 
(source 15 in lSwartz et alj [2003). However, the 21cm radio source 
is too large for an SNR of the required age and there is no convinc- 
ing supporting evidence of a source within the other images stud- 
ied. Using the published known X-ray sources we might have ex- 
pected to find an SNR if it was very young or bright and we would 
have expected to find an associated SNR using the Ha images. We 
identified no convincing associated SNR candidates within the er- 
ror ellipse. 

Although it has been determined that the error ellipse does 
cross potential star formation regions as required by the majority of 
SGR models, it should also be noted that this is not essential for all. 
An alternative route has been p roposed for producing a magnetar 
by white dwarf (WD) me rgers dKing. Pringle, & Wickramasinghel 
1200 ll ; iLevan et alj|2006bh . As WD have long lifetimes, WD-WD 
mergers would be associated with older populations of stars. It is 
possible that accretion induced collapse (AIC) will drive off a frac- 
tion of the envelo pe, leaving something akin to an SNR behind (e.g. 
iBaron et al.ll987t) . The mechanisms underlying AIC are poorly un- 
derstood, and the physical characteristics and detectability of such 
remnants are not clear. Therefore a SGR produced through this 
channel could be formed in an old stellar population within the 
outer disk or halo, and the non-detection of an SNR within the re- 
gion does not place constraints on this model. 

If the progenitor was an SGR giant flare, then there might 
be significant similarities in the light curve and spectrum of GRB 
051103 to the giant flare from SGR 1806-20. lofek et alj d2006h 
noted that the light curve of these two events were consistent, i.e. 
the light curve of GRB 051103 is similar to what would be ex- 




Figure 6. This shows an image taken using GALEX showing the UV sources 
within the error ellipse. Two of the brightest sources discussed bv lOfek et all 
12006), within the ellipse, are highlighted by the red circles. 



pected from an extragalactic version of the giant flare from SGR 
1806-20. In Table 1, we have shown, for the joint RHESSI + KW 
fits, that initially a = ^ and softens to a = 0.39+° jo. 

Although this is unusually hard for a GRB, it is consistent with the 
photon index of the giant flare from SGR 1806-20, a = 0.2 ± 0.3 
dPalmer et aT]|2005h . The peak luminosity of GRB 05 1 103, assum- 
ing it was from an SGR in M81, is approximately 4.7 x 10 48 erg 



s . This is a factor of 10 brighter than the peak luminosity of 
the giant flare from SGR 1806-20, w hich is 2 - 5 x 10 47 erg s" 1 



assuming it is at a distance of 15kpc i Hui'lev et al. 2005). With the 
revised distance estimate from Bibby et al. (2008), the peak lumi- 
nosity of the giant flare from SGR 1806-20 would be 7 x 10 46 erg s 
~ , suggesting that a much smaller percentage of SGRBs are SGR 
giant flares. This value is 30 times fainter than the peak luminos- 
ity of GRB 051103 if it was from an SGR giant flare in M81 and 
in this case GRB 051103 would be the most luminous SGR giant 
flare observed. In comparison, the peak luminosity of GRB 070201 
is 1.14 x 10 47 erg s _1 assuming it was in M31 dOfek et alj|2008h . 
which is an order of magnitude fainter than GRB 051 103 and com- 
parable to the giant flare from SGR 1806-20. It is important to note 
however, that there is c urrently no theoretical upper limit for the 
energy of a giant flare. iDuncan & Thompson! dl992h showed that 
the total energy available is given by E oc 3 x 10 47 B 2 5 erg 
where B15 = B/10 15 G. Therefore, the magnetic dipole (B) of 
SGR 1806-20 would only need to increase by a factor of ~5 to 
produce a giant flare with an energy that is 30 times greater than 
the one from SGR 1806-20. 

Although the gamma-ray data suggest that GRB 051103 may 
be an extragalactic SGR giant flare, it is important to note that 
SGR giant fla res are rare events. Consider i ng pl ausible luminos- 
ity functions, IChapman. Priddev. & Tanvirl d2009i) calculated the 
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probability that the IPN would observe a giant flare, with energy 
greater than the energy emitted by the giant flare from SGR 1806- 
20, in the region surveyed during the 17 years it has operated. For 
one giant flare, they calculated the probability to be 10%. How- 
ever, as we discussed in the introduction, there are four poten- 
tial candidates for extragalactic SGR giant flares, including GRB 
07 0201 near M31 whic h has been identified as an SGR giant flare 
bv lMazets et al.l J2008h . The probability that the IPN has detected 
two SGR giant flares, wit h energy greater than the giant flare from 
SGR 1806-20, is 0.6% dChapman, Priddev, & Tanvhl 120091) . Re- 
cently, several new SGR candidates have been ide ntified including 
0501+4516, 1550-5418 and possibly 0623-0006 jBarthelmv et all 
l2008ahlKrimm et alj20 08; Barth elmv etal]2 008b), which may im- 
ply that the number of SGRs in the Milky Way is higher than pre- 
viously thought. In this case, the luminosity of the giant flare from 
SGR 1806-20 would have to be at the peak of the luminosity func- 
tion of SGR giant flares and therefore giant flares of this luminosity 
must be extremely rare events. This argues that GRB 051 103 is un- 
likely to be a second SGR giant flare in the nearby universe. 



move the requirement for an SNR. Additionally, the luminosity of 
GRB 051103, assuming it is from an SGR giant flare in M81, is 
significantly higher than known SGR giant flares but still attainable 
with current theoretical models. Giant flares with luminosity simi- 
lar to the giant flare from SGR 1806-20 are extremely rare and it is 
unlikely that GRB 051103 and GRB 070201 are both extragalactic 
SGR giant flares. 

Although we have not considered this option in detail, it is 
possible that the progenitor of GRB 051103 was a compact binary 
merger in M81. In this case, it would just be within the reach of 
current gravitational wave searches. This scenario was ruled out 
at >99% confidence for GRB 070201 in M31 using the Laser In- 
terferometer Gravitational- wave Observatory (LIGO) observations, 
and distances out t o 3.5 Mpc were ruled out to 90% confidence 
dAbbott et all 120081) . The LIGO Scientific Collaboration is cur- 
rently considering a search for gravitational-wave signals in the 
data surrounding GRB 051103 (G. Jones and P. Sutton, private 
communication). 



6 CONCLUSIONS 

GRB 051103 illustrates the difficulties of identifying a short burst 
as an extragalactic giant magnetar flare beyond a reasonable doubt. 
Even setting aside the questions of detecting and localizing such 
events, and establishing their associations with nearby galaxies, 
their interpretation is problematic. On the one hand, the localiza- 
tion, short duration, and hard energy spectra of GRB 051103 sug- 
gest that it is a giant flare from M81. However, a deeper analysis of 
its time history demonstrates that the periodic component, which 
is a key signature of giant flares, is unlikely to ever be detected at 
great distances by the IPN if all giant flares are similar to the three 
observed to date. The energy spectrum at the peak of the emis- 
sion is very hard (E pea k ~3 MeV), and is detected to 7 MeV at 
the 3<t level, with marginal emission up to 17 MeV. Yet it is not 
inconceivable that a short duration GRB could have these proper- 
ties. Although the E pea fc of GRB 05 1 103 evolves from hard to soft, 
the evolution to a very soft spectrum, which is expected during the 
oscillatory phase of an SGR giant flare, is undetectable, as is the 
oscillatory phase itself. Thus evidence for an extragalactic giant 
flare origin of GRB 051103 remains tantalizing, but inconclusive. 
On a more positive note, if an extragalactic magnetar flare occured 
within the Swift-BAT field of view, so that the XRT could begin 
observing within a minute or so, the periodic component would be 
detectable at low energies to at least 10 Mpc (Hurley et al. 2005). 

We have presented new optical observations of GRB 051103 
and have determined that there is no R band optical afterglow with 
a limiting magnitude of ~ 24.5 (for an afterglow overlapping a host 
galaxy, the limiting magnitude is ~ 23.3) in the region of the er- 
ror ellipse covered by our observations. Comparison of the prompt 
emission of GRB 05 1 103 with a sample of other SGRBs leads us to 
conclude that if it was a classical SGRB we would expect to have 
located an optical afterglow in our observations. 

In contrast, if GRB 051103 were an SGR giant flare in M81, 
non-detection of an afterglow would not be surprising as the ex- 
pectations for optical afterglow emission lie significantly below the 
limits obtained here, or the limits likely to be attained via current 
technology. The case for an SGR origin would be strengthened 
if there were an accompanying SNR within the error ellipse, but 
there is no evidence of t his. An SGR produ ced via accretion in- 
duced collapse of a WD dLevan et alj|2006bh would, however, re- 
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